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Abstract 

Quantum chaos — the study of quantized nonintegrable Hamilto- 
nian systems — is an extremely well-developed and sophisticated field. 
By contrast, very little work has been done in looking at quantum 
versions of systems which classically exhibit dissipative chaos. Using 
the decoherence formalism of Gell-Mann and Hartle, I find a quantum 
mechanical analog of one such system, the forced damped Duffing os- 
cillator. I demonstrate the classical limit of the system, and discuss 
its decoherent histories. I show that using decoherent histories, one 
can define not only the quantum map of an entire density operator, 
but can find an analog to the Poincare map of the individual trajec- 
tory. Finally, I argue the usefulness of this model as an example of 
quantum dissipative chaos, as well as of a practical application of the 
decoherence formalism to an interesting problem. 



1 Introduction 
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1.1 Classical laws and dissipative chaos. 

Recently, Gell-Mann and Hartle, among others [Jj, 0, |!|, have studied the 
problem of classical laws arising from quantum theory in the light of the 
decoherence formalism. In this approach, one considers possible histories of 
a given system to which probabilities can be assigned that obey classical 
probability sum rules. In order for histories to decohere in this way, it is 
usually necessary to coarse-grain the description of the system, by giving the 
values of its variables only at certain times, or averaged over certain intervals, 
or by neglecting certain variables and retaining others, or a combination of 
all of these. 

They have found that within this formalism it is possible to define in a 
very rigorous way the classical equation of motion based only on the underly- 
ing quantum theory. In doing so, both dissipation and noise typically appear, 
arising as a consequence of coarse-graining over neglected degrees of freedom. 
This takes advantage of the well-studied phenomenon of environmentally in- 
duced decoherence || . In addition to casting new light on the problem of how 
classical laws of physics arise, this provides an unparalleled tool for study- 
ing quantum mechanical systems with dissipation, and seeing how this alters 
their behavior from the more usual Hamiltonian behavior of closed systems. 
One area which can profitably be treated this way is dissipative chaos. 

There has been an enormous amount of work done on "Quantum Chaos," 
i.e., quantizing nonintegrable Hamiltonians which classically exhibit chaotic 
behavior. This has turned up beautiful connections between classical chaotic 
behavior and their quantum quasiperiodic equivalents. But very little has 
been done in looking at the quantum versions of systems which classically 
exhibit dissipative chaos, or on looking at their classical limit [TJJ. Classical 
dissipative chaos is qualitatively very different from Hamiltonian chaos, and 
one would expect their quantum equivalents to reflect this difference, but this 
has not been widely investigated. Indeed, even very extensive treatments of 
quantum chaos rarely deal with dissipative systems at all 0. 

There are a number of reasons for this. The first is that dissipation is diffi- 
cult to treat in normal quantum mechanics. The usual Schrodinger equation 
is only valid for closed systems without friction. Open systems in general, 
and dissipation in particular, can be handled using the influence functional 
approach of Feynman and Vernon || ; this has been done in the case of Brow- 
nian motion by Caldeira and Leggett J7| among others. This approach has 
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not been widely used, though, until recently, as it involves considerable con- 
ceptual and mathematical baggage || . Also, the types of behavior of most 
interest to those who study quantum Hamiltonian chaos involve the coher- 
ent evolution of the wave function, with its attendant complicated structure 
(e.g., the "scarring" of energy eigenfunctions about classical periodic orbits, 
the statistics of energy level spacing). The presence of strong damping wipes 
out this coherent structure. 

Chirikov et al. have summed up the usual attitude towards quantum dis- 
sipative chaos: "In what follows we will discuss only Hamiltonian (nondis- 
sipative) systems, considering them to be the more fundamental ones. Phe- 
nomenological friction is but a crude approximation of the molecular Hamil- 
tonian chaos which is inevitably related to some noise according to the 
fluctuation-dissipation theorem." And further, they divide the problem of 
quantum chaos into two parts, the quantum dynamics of the wave function 
in isolation, and the results of measurement "with its unavoidable statistical 
effect of the irreversible ip collapse which is a sort of inevitable noise." || 

While this is undeniably true, most systems are not isolated, and so it is 
perhaps useful to consider systems for which dissipation is important. Dis- 
sipative chaotic systems may not be fundamental, but they are nevertheless 
interesting. Decoherence is an appropriate formalism in which to study them 

0- 

In the rest of this section I give a brief introduction to the decoherent 
histories formalism of Gell-Mann and Hartle. Then in section 2 I derive 
a model for a quantum forced, damped nonlinear oscillator, following the 
usual system/environment coarse- graining. In section 3 I discuss the classical 
properties of the forced, damped Duffing oscillator, and describe some of the 
properties of dissipative chaos which it exhibits. In section 4 I treat the 
quantum version of this problem, and show how one can make close contact 
with the classical theory using the decoherent histories formalism. In section 
5 I illustrate this with a numerical example, and in section 6 I summarize 
my conclusions. 

1.2 Decoherent histories. 

In the formalism of decoherent histories, systems are described by a set of 
exclusive and exhaustive histories {a}, which can be thought of as different 
possibilities for the system's evolution. While there is a vast range of possible 
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sets of histories to choose from, these sets are restricted by the decoherence 
condition 

D[a,a'] = p a 5 aa ,, (1) 

where D[a, a'] is the decoherence functional and p a is the probability of 
history a occurring. This decoherence condition restricts one to histories 
which obey the usual probability sum rules; these histories do not interfere 
with each other. 

To make this more explicit, in ordinary nonrelativistic quantum mechan- 
ics one can specify a history by enumerating a complete set of orthogonal 
projection operators {PL(ti}} at a sequence of times U. A single history 
is then given by choosing one projection operator at each time. This is 
equivalent to enumerating a set of possible assertions about the system at 
a sequence of times, and having each history be a string of such assertions. 
One can define a history operator 

C a = PZ(t n )...Pl(t 2 )P 1 ai (t 1 ), (2) 

where a is a shorthand for the choices at times tj. The decoherence 
functional is then 

D[a,a'] = Tr{C aP Cl,}. (3) 

The density operator p is the system's initial condition. 

As a rule, it is impossible for very fine-grained histories to decohere; 
thus, considerable coarse-graining is required. One very common coarse- 
graining used to study decoherence in systems with many degrees of freedom 
is to completely trace out certain freedoms (the "environment" ) while leaving 
others completely fine-grained (the "disinguished subsystem" ) . This was first 
studied by Feynman and Vernon || and applied to decoherence by Zurek || 
among others. We will initially be considering this type of coarse-graining. 

2 The Model 

The particular model we will study is based on earlier work on decoherence 
in systems with dissipation || 0], @, I n this model we will divide our 
system into a distinguished variable x, termed the system variable, and a 
set of reservoir variables {Qk} which we will trace over. This system and 
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reservoir will have a total action 

S[x(t), Q(t)] = S sys [x(t)} + S ICS [Q(t)} + f M„t(x(t), Q(t))dt, (4) 

J to 

where the system variable will be treated as a particle moving in a potential 

S sys [x(t)} = jf" (y* 2 (0 - U(x(t))yt, (5) 

the reservoir is approximated as a collection of harmonic oscillators 

sUQ(t)] ^E/^ftW 2 -*^) 2 )^ (6) 

and the interaction is linear in x and Q: 

Vmt (a;, Q) = -x IkQk- (7) 

k 

We will make the additional assumption that the initial density matrix of 
the system and reservoir factors, and that the reservoir is initially in a ther- 
mal state. Then p t otai(^> Q',x', Q') = x( x l ^OV'oCQj Q')> where i/; — Pt is a 
thermal density operator at temperature T. 

The decoherence functional in this coarse-graining is then 

D[x'{t),x{t)\ =exp^-! [ S sys [x'(t)} -S sys [x{t)\ +W[x'{t),x{t)}\ X {xo;x' ). (8) 

W[x'(t), x(t)} is the influence phase, which includes the collective effects of 
the traced-over reservoir degrees of freedom. As shown by Caldeira and 
Leggett 0, this functional is 

■ 2 

W[x'(t), x(t)\ = y coth( WfcT) 

x dt ds cos(u k (t — s))(x (t) — x(t))(x'(s) — x(s)) 

Jto Jto 

^— f * dt I ds sin(cj k (t - s))(x'(t) - x(t))(x'(s) + x(s)). (9) 

2mUk Jto Jt 
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We can now switch to new variables X = (x + x')/2 and £ = x — x'. In 
these variables 

S sys [x'(t)]-S sys [x(t)) = 

£' (-Mi(t) - ^(x(t)j)dt - £omx + o(e), m 

which as we see contains the Euler-Lagrange equations. We also go to a con- 
tinuum of oscillator frequencies with a Debye distribution, in which the dis- 
crete sums become integrals over a weighting function g(u>) = rju 2 exp(— u/fl), 
where Q is a large cutoff frequency, so that 1/Q « tf — 1 . In this limit, the 
influence phase becomes 

w[x(t),m = £ / £(t)(-MAX(t)-2rx(t))rft 

+iK f S e(t)dt + 0(n- 2 ). (11) 

J to 

where A = r]Q/m, T = uri/AmM, and K = 4MTkT/h. The A term has the 
form of a linear force; it can be absorbed into the system action by going to 
an effective potential 

U cS (X) = U(X) + MAX 2 /2. (12) 

The T term has the form of a dissipation. 

The imaginary term is of particular interest. It suppresses D[X(t),£(t)] 
when £ ^ 0. Since £ ^ corresponds to the "off-diagonal" terms of the 
decoherence functional (x(t) ^ x'(t)), the suppression of these terms results 
in approximate decoherence of this set of histories. This suppression of off- 
diagonal terms is clearly related to the presence of noise || |]]. The kernel 
of this term can be identified with the two-time correlation function of a 
stochastic driving force F(t) in the quasiclassical limit. This correlation 
function is 

(F(t)F(s)) = %K8(t - s) (13) 

in the continuum case, with (F(t)) = 0. So in the quasiclassical limit this 
system obeys the classical equation of motion 

t + ^^W + ZT^m/M. (14) 
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Instead of taking the reservoir to be in a thermal state initially, we can 
take it to be in a displaced thermal state, 

Pdt = D(q(uj),p(uj))p T b(q(uj),p(uj))\ (15) 

D(q(u),p(u)) is the coherent state displacement operator: 

D(q(uj),p(uj))\0) = \q(uj),p(uj)), 

where \q(u),p(w)) is the coherent state centered on (q(uj) , p(u>)) at frequency 
uj. If we take q{ui) to be sharply peaked around a certain frequency, q{uj) = 
q5(u — uj ), and p(u) = 0, then in the quasiclassical limit the above equation 
of motion QUI) gains an additional term 

x + —^L(x) + 2Tx = qcos(tu t) + F(t)/M. (16) 
M ax 

This is exactly the form for a nonlinear oscillator with damping and a periodic 
driving force, with additional noise. 

In a truly classical system, F(t) would vanish as T — > 0, but in the 
quantum theory noise is always present, even at absolute zero. One can think 
of it as arising from the zero-point oscillations of the reservoir oscillators. At 
low temperatures, however, the two-time correlation function of the noise is 
highly nonlocal in time. At T = 0, 

Re W[X(t),£(t)] ~ f tf dt f f ds £(t)t(8)/(t - s) 2 . (17) 

Correlations in the noise persist for all times. Because of this form of the 
kernel, doing exact (or even numerical) calculations in the low-temperature 
limit is extremely difficult. This is why the high T limit is generally used. 

3 The Classical Forced Damped Duffing Os- 
cillator 

We are interested in finding quantum equivalents to classical systems which 
exhibit dissipative chaos. While many such systems (e.g., fluid mechanics) 
have no easily realizable quantum limit, there are some which can be readily 
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quantized as shown in section 2. These are the nonlinear oscillators with 
damping and driving. 

One much-studied classical nonlinear oscillator is the Duffing oscillator, 
characterized by a two-welled polynomial potential, 

4 2 

if* T* 

U{x) = X - - X -. (18) 
With forcing and damping, this gives an equation of motion 

x + (x — x 3 ) + 2Tx = q cos(ujot), (19) 
where we take M = 1. This system is chaotic for certain values of q, T, and 



Uq [O. For example, q = 0.3, T = 0.125 and co>o = 1.0 is a common choice. 

Since this system has explicit time dependence, its phase space is three- 
dimensional, (x,p,t). It is common to discretize the dynamics by taking a 
Poincare section, considering only the points on a surface of constant phase 
(xi,pi) at times £, = lixijuj^. The continuous dynamics defines a map f: 

(xi,pi) -> (zi+ijPi+i) = f(:Zi,Pi). (20) 

If the oscillator is non-chaotic, there is a stable attracting fixed point or group 
of periodic points to which the (xi,pi) quickly tend. These points correspond 
to a periodic orbit of the continuous dynamics. When the oscillator becomes 
chaotic, the stable set becomes a strange attractor, a fractal structure with 
non-periodic behavior. There are, in addition, an infinite number of unstable 
fixed points and periodic points. (See figure 1.) 

We can also look at the classical dynamics from the point of view of 
probability measures P(x,p) on phase space. The map f of points in phase 
space induces a map on probability measures 

- P, +1 (x,p) = *((,,„) - f(*V))P,(*V). (2D 

By means of this sort of map we can readily make contact with the quantum 
theory. 

Of particular interest are invariant probability measures Pj nv . There are 
many of these, most corresponding to unstable fixed points and periodic 
points of the map f. It is possible to eliminate these unstable solutions by 
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including a small amount of noise in the equation of motion (|T^) . This effec- 
tively broadens the delta functions in (21) into peaks of finite width e, and 
eliminates all unstable solutions, leaving a single unique P inv corresponding 
to the strange attractor. Classically, we can then allow the noise to go to 
zero, and look at Pi nv in the zero-noise limit. In that limit, the invariant 
probability measure becomes a generalized function with substructure at all 
length scales. It is a fractal. 



4 Decoherent Histories, Quantum Maps, and 
Probability 

At best, the functional described in (|8|) can only be approximately decoher- 
ent. Clearly, off-diagonal terms will not vanish for sufficiently small 
More coarse-graining is needed in the description of x(t) and x'(t). Also, 
specifying a value, even an approximate value, of x(t) for all times t is an 
extreme fine graining. It is more common to instead specify x at a series 
of discrete times U. Thus, instead of a complete trajectory x(t) one gives 
only a series of x values {x{\. Coarse-graining in position as well, one could 
divide up the range of x into finite non-zero intervals A*- , where j is an index 
specifying which interval x fell in at time £j. A history would now be a se- 
ries of indices {«;}, specifying that x fell in the interval A^. at time £j. Note 
that to achieve decoherence, these times U cannot be too close together; they 
must generally be separated by at least the decoherence time M. For high 
temperature systems this is typically quite short, of the order % /2MTkTd 2 , 
where d is the size of the intervals. 

Such a coarse-graining gives us a new decoherence functional: 

D[a, a') = / 5x [ 5x' D[x(t),x'(t)\, (22) 

J a J a' 

where the limits specify integration only over those paths which pass through 
the series of intervals A^. at the times tj. The probability of a given history a 
is of course given by the diagonal terms of this functional. Since the original 
decoherence functional given by (|8|) has an exponent quadratic in £, the path 
integrals over £ can be carried out; we then let a = a' and get 

p{a) = Jj^JjX exp|--^-£e 2 (t)rft}«;(Xo,MAo), (23) 
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where e(t) = MX + (dU cS /dt){X) + 2MTX - qcos{u t) is the right-hand 
side of the equation of motion |10j . From this we see that the probability will 
by peaked about histories which approximately obey the classical equation 
of motion e(t) = 0, more and more sharply as we approach the classical limit 
where M is large. The w(X , MX ) is the initial Wigner distribution of the 
system. 

The Wigner distribution is defined in terms of the density matrix: 

w(X,p) = - f e-^ h X (X + e/2; X - £/2) df. (24) 

7T J 

The distribution behaves very similarly to a classical probability distribution 
on phase space, except that w(X,p) can be locally negative (though it must 
sum to 1 over all of phase space, and be non-negative on average over regions 
with volumes larger than h). The expectation values of functions of X and 
p can be calculated by averaging them over phase space using w(X,p) as a 
weighting function, though there is usually some ambiguity about the order- 
ing of operators. As one goes to the classical limit, on scales large compared 
to h, this ambiguity becomes unimportant. 

An interesting way of looking at this system is in terms of the evolution of 
the Wigner distribution with time. If we consider surfaces of constant phase, 
as in the classical case, we can define a quantum map, 

Wi — > w i+ i = Twi, 

u>;+i(Xi,pi) = J dX J dp T(X 1 ,p 1 ;X ,p )w(X ,p ). (25) 
The transition matrix T is defined by the path integral 

T(X 1 ,p 1 ;X ,p ) = - f dCod^e-^-^ 

7T J 

x J 6X6£ exp^S sys [X(t)+C(t)/2] 

-S sys [X(t) - e(t)/2] + W[X(t),t(t))}, ( 26 ) 



^o^ie~ i(Cm " ?oPo) 
1 rU+i 



xj'5X exp j-— jf 1+1 e 2 (t)dt + i^M^ - ^ MX )}, (27) 
5X 5(po - MX )5( Pl - MXi) exp[-— J e 2 (t)dtj. 
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This evolution strongly resembles the classical evolution of probability mea- 
sures induced by the phase-space map, and in the classical limit we expect 
w(X,p) to evolve towards an invariant distribution w inv (X, p) which closely 
resembles the classical invariant measure Pi rcv (x,p). Graham |4| has demon- 
strated this sort of behavior in his work on the quantum Lorenz model, which, 
though very different in approach from this paper, may nevertheless be in- 
dicative; and the author's own numerical simulations |H| seem to bear this 
out (though of course one would not expect numerical simulations to exhibit 
unstable alternative solutions). 

In the quantum case, it is impossible for the noise to ever truly vanish. 
Even at absolute zero, zero-point fluctuations remain that prevent w(X,p) 
from becoming a true fractal. Though the invariant distribution may strongly 
resemble P imr for a wide range of scales, there is always some scale at which 
the quantum noise "smears out" w inv (X,p). 

While these maps on Wigner distributions make contact with the classi- 
cal theory, ideally we would like to find some quantum analog to (^3), i.e., 
a description in terms of individual histories, rather than probability distri- 
butions. To do this, let us consider yet another coarse graining. Consider 
the decoherence functional (^), where we take the sequence of times to be 
those corresponding to the surface of section t« = 2ni, and let the projection 
operators P QtP be onto localized cells of phase space centered at (q,p). While 
there are no true projections onto cells of phase space, there are approximate 



projectors which can be used to get approximate decoherence |T2|, |10| . For 
example, simple coherent state projections \q,p){q,p\ can be used. 

A history can then be specified by a series of points {qi,Pi} at the times 
ti, and the decoherence functional calculated 

D[{qi,Pi}, WM] = Tr{P qntPn T(. ■ ■ T(P guPlPo P q[>p[ ■ ■ -)P gk , pli }. (28) 

Here we have taken T to be the transition matrix on density operators rather 
than Wigner distributions; it is simple to go from one representation to the 
other. This is a quantum surface of section. At each time ti the system is 
localized in a cell in phase space centered on (qi,Pi), and probabilities can 
be assigned to each possible next point (q i+ i,p i+ i). This differs, of course, 
from the classical case where the evolution is deterministic; but from ( p3|) it 
is clear that these histories will be peaked about the classical evolution in 
the quasiclassical limit. This is shown explicitly by the numerical example 
in the next section. 
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5 Numerical simulation and quantum state 
diffusion 



While the decoherent histories formalism has great interpretational power, 
it is not very convenient for numerical simulation. Enumerating all the pos- 
sible histories and calculating the elements of the decoherence functional is 
a daunting task. What one would like is a method of generating the his- 
tories with the correct probabilities without having to solve the full master 
equation at every step. 

Recently, it has been shown that the theory of quantum state diffusion 
provides just such a technique. Quantum state diffusion is one of the so- 
called quantum trajectory methods, in which the master equation evolution 
is unravelled into quantum trajectories of individual states. These states obey 
a nonlinear stochastic differential equation which in the mean reproduces the 
master equation. Because one need deal only with a single quantum state at 



a time, this is very suitable for numerical calculation [16 



Diosi, Gisin, Halliwell and Percival [17] have shown that these individual 
quantum trajectories correspond to a set of approximately decoherent his- 
tories. In the case of a dissipative interaction, these correspond to histories 
of systems localized into small cells in phase space, and their probabilities 
match those given by decoherent histories. Thus it is ideal for the sort of 
problem we are interested in. 

For further details see the references. More work on the connections be- 
tween decoherent histories and quantum state diffusion, and their application 
to dissipative chaos, is currently underway |13| . 

In figure 2 we see one such trajectory, generated in the quasiclassical limit 
(where h = 10~ 4 ). One can see that this is very close to the classical limit, 
but with additional "smearing" due to the presence of noise. This smearing 
sets a lower cutoff scale to the substructure of the strange attractor. As one 
continues to go to the classical limit, more and more substructure appears, 
and the noise becomes less and less important. 

Note that in this chaotic system we expect the trajectory to evenly sample 
the "invariant" Wigner distribution W{ nv over time. From the distribution of 
points we see that this is indeed very close to the structure of the classical 
strange attractor in figure 1. 
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6 Conclusions 



As we have seen, it is possible to use the decoherence formalism to study at 
least some systems exhibiting classically chaotic behavior, and to do so in a 
way which includes dissipation in a simple and natural fashion. Though the 
model treated here is very much a special case, intended only to illustrate the 
basic ideas of the theory, it is remarkable how many details can be brought 
out and studied by its means. 

Certainly these techniques should work for any kind of nonlinear oscil- 
lator, or for multivariable extensions of them. It might well be possible to 
treat systems of experimental interest, arising in fields such as quantum op- 
tics. Some work on such systems has already been done by other workers 

B 

Using the usual master equation formalism it is possible to draw a close 
connection between the classical theory of probability measures and the quan- 
tum Wigner distribution. But with decoherent histories, one can also find a 
quantum analog to individual chaotic orbits, such as the quantum surface of 
section defined in section 4. 

One can then argue analytically that these quantum histories become 
more and more closely peaked about the classical equations of motion as one 
goes to the classical limit; and this correspondence can also be demonstrated 
numerically. 

Further analytical study may yield better results for the probabilities 
and decoherence of phase space histories. And it may be fruitful to explore 
what equivalents there are in the quantum case to classical quantities such 
as Lyapunov exponents, fractal dimension, and Kolmogorov entropy. This 
theory should amply reward further study, both analytical and numerical. 
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Figure 1. The classical forced damped Duffing oscillator surface of section 
in the chaotic regime, q = 0.3, T = 0.125, uq = 1.0. 



Figure 2. The quantum forced damped Duffing oscillator surface of sec- 
tion, generated by the quantum state diffusion algorithm in the quasiclassical 
limit, k = 10" 4 , q = 0.3, T = 0.125, cu = 1.0. 
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